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Abstract 

We have made a new calculation of the flux of secondary positrons above 
100 MeV expected for various propagation models. The models Investigated are 
the leaky box or homogeneous model, a disk-halo diffusion model, a dynamical 
halo model and the closed galaxy model. The parameters of these models have, 
In each case, been adjusted for agreement with the observed secondary/primary 
ratios and ^%e abundance. The positron flux predicted for these models Is 
compared with the available data. The possibility of a primary positron 
component Is considered. 
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INTRODUCTION 

Low energy positrons are expected to be produced by the decay of 
radioactive Isotopes created by nucleosynthesis In supernovae (e.q., Colgate, 

1970) and possibly by pair production near the surface of pulsars (Sturrock, 

1971) . The bulk of the cosmic ray positrons observed above 100 MeV are, 
however, thought to be of secondary origin, resulting from the decay of 
produced In nuclear interactions of cosmic rays In the Interstellar medium. 
Observations of high energy cosmic ray positrons, when combined with model 
predictions may thus help us understand the propagation of cosmic rays In the 
galaxy. Previous calculations of the production rate of positrons from this 
source have differed by as much as 50 percent resulting In diverse conclusions 
regarding propagation. In the present paper, we give the results of a new 
calculation of the positron production rate and estimate the flux of cosmic 
ray positrons expected for various propagation models. 

MODELS FOR COSMIC RAY PROPAGATION 


The propagation models we will consider include the widely used leaky-box 
model, a conventional diffusion model, the dynamical halo model and the closed 
galaxy model. In each case the model parameters are adjusted so that the 
predicted energy dependence of the boron/carbon ratio and the surviving 
fraction of the radioactive nuclide i°Re are consistent with observation. 

Other propagation models will be briefly discussed. 

a) Leaky Box or Homogeneous Models 


These models are characterized by an exponential distribution of 
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cosmic ray ages with mean lifetime, <t>, which is related to the mean escape 
length or 'grammage', Xg, by 

<t> = Xg/pBC (1) 

where p is the mean density of interstellar material as sampled by the cosmic 
rays. From boron/carbon and other secondary to primary ratios, Protheroe, 
Ormes and Comstock (1981) found that Xg « 7 g/cm^ at rioidities, R, less than 
4 GV/c and Xg « g/cm^ above 4 GV/c. Recent satellite 

measurements of the abundances of isotope, of Be by Wiedenbeck and Greiner 
(1980) indicate that ~ 29 percent of l^Be survives decay at interstellar 
energies of a few hundred MeV/nucleon. From this surviving fraction, 
Wiedenbeck and Greiner conclude that <t> « (B^4^2*4 ^ ^ years for relati- 
vistic particles. For a mean escape length of 7 g/cm^ of interstellar matter 
(90 percent hydrogen and 10 percent helium by number) this corresponds to a 
mean density of interstellar nuclei, n, of (0.40 ) atoms/cm^. 

b) Diffusive Halo Model 

In this model, cosmic ray sources and matter are located in a 
disk of thickness 2a, which is surrounded by a halo of thickness 20. Cosmic 
rays diffuse throughout the disk and the halo and escape freely from the 
boundary of the halo (Prischep and Ptuskin, 1975). We adopt a one dimensional 
approximation similar to that of Owens and Jokipii (1977a), and assume that 
the diffusion coefficient in the halo has the same value as that in the disk. 
From Freedman et a1. (1980), we find that the diffusion coefficient, k, is 
proportional to 3/Xe» constant of proportional ity depending on the values 
of a, 0 and the matter density in the disk. We adopt a « 100 pc, comparable 
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to the half thickness of the total hydrooen* 1.e., n(HI) + ?!n{H2K layer and a 
density of 1.1 atoms/cin^ correspondino to a hydrogen surface density of ^.0 
M0/pc^p the local value (Gordon and Burton, 1976). For an interstelar mediuin 
containing 10 percent He by number, this corresponds to a surface density of 
7.2 Mq/PC^. The halo is assumed to be devoid of matter. Since Ag « 7 q/cm^ 
at low energies, we can Obtain values of the diffusion coefficent, k, which 
would result in this grammage for different values of the halo thickness, n. 
This relationship may then be combined with the observed surviving fraction of 
^°Be, fj, to fix both D and < in this model. Using formulae of Freedman et 
al. (1980) and fg** 0.29 ± 0.Q8 (Wiedenbeck and Greiner, 1980), we obtain 0 = 
(1.7^Q*g) kpc and k/b = (I.G^q*^) x 1020 cm2/s (note that the errors in n and 
ic are not independent). This halo size is consistent with that obtained from 
an analysis of gamma ray data by Stecker and Jones (1977). Above a rigidity 
of 4 GV/c we adopt k/3 = 1.6 x 1020 (r/4)0.4 cm2/s. 

c) Dynamical Halo Model 

This model is similar to the diffusive halo model described 
previously except that cosmic rays are convected outward in the halo by a 
galactic wind (Jokipii, 1976). The velocity of the scattering centers, or 
convection velocity, is assumed to be zero in the disk and a constant value, 

V, in the halo. Again, we assume that the value of the diffusion coefficient 
in the halo is the sanne as that in the disk. Because of the outward 
convection, a larger halo is required in this model to fit the observed values 
of Ag and fg than for a static halo. The surviving fraction of ^^Re in the 
dynamical halo model has been discussed by Owens and Jokipii (1977a), Jones 
(1979) and Freedman et al. (1980). Particles diffusing across the disk-halo 
boundary lose energy by shock deceleration and so Ag and fg depend on their 
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energy spectrum In addition to the parameters of the propagation model. The 
motivation for this model comes from the observed energy dependence of the 
grammage. In this model, the diffusion coefficient may have a power law 
dependence on rigidity at all energies giving the observed decrease vdth 
energy above a few GeV/nucleon, but still give Xg - constant at low energies 
because of the galactic wind. Jones (1979) found that the form, 

ic = SKq R^-5 (2) 

gave a good fit to the observed energy dependence of Xg at high energies. 

At low energies, if the cosmic ray injection spectrum of secondary nuclei 
had a differential power law exponent, of 2.5, the observed energy dependence 
of Xg was also well fitted by this model for VD/kq “ (GV/cf'^^and pac/V = 

20 g/cm^. The observed spectrum at low energies is however more consistent 
with an injection spectrum of the form: (T + 400 MeV/nucleon)"^*^, where T is 
the kinetic energy per nucleon (Garcia-Munoz et al., 1977) rather than with 
T"^*^. At a few hundred MeV/nucleon this spectrum is similar to the form f?T”T 
for Y ~ 1.65. Using the results of Freedman et al. (1980) for this value of 
Y, we find that better agreement with Xg at low energies may be obtained with 
pac/V « 14.4 g/cm^, or V « 15.3 km/s. To obtain D (and Kq) we again use the 
observed surviving fraction of ^°Be. The formulae of Freedman et al. (1980) 
yield: D = (4.0 kpc; k q= (l.^^J;^) x lO^V/s. 

d) Closed Galaxy Model 

Electrons and positrons have been considered in the closed 
galaxy model of Rasmussen and Peters (1975) by Radhwar and Stephens (1976), 
Ramaty and Westergaard (1976) and French and Osborne (1976). Problems with 


this model led to Its revision by Peters and Westergaard (1977) and this Is 
the model we shall consider here. The Inability of conventional propagation 
models to explain the high cosmic ray antiproton flux observed at ~ 10 GeV 
(Golden et al.* 1979) has led to a resurgence of Interest In this model 
(Protheroe, 1981; Stephens, 1981). Secondary positrons In this model have 
been considered in an approximate way by Gller, Wdowczyk and Wolfendale (1977) 
and Stephens (1981). 

In the closed galaxy model of Peters and Westergaard (1977) cosmic ray 
sources are located In the spiral arms of the galaxy. Cosmic rays are then 
partially trapped In the arms and leak out slowly Into the surrounding halo, 
the outer boundary of which constitutes a closed box from which they cannot 
escape. Depletion of cosmic ray nuclei in the halo which contains low density 
Interstellar matter is then due solely to nuclear interactions and energy 
losses. The halo thus contains an "old component" of cosmic rays consisting 
mainly of protons; heavier nuclei leaking from the arms spall into nucleons. 

In this model the Sun is located in a spiral arm and the cosmic rays we 
observe comprise a "young component" from the sources (these cosmic rays have 
not yet escaped from the arms) plus the old component which permeates the 
whole galaxy. 

The parameters describing the closed galaxy model are K, the ratio of the 
mass of interstellar material in the galaxy as a whole to that in the spiral 
arms, and nj^, the number density of interstellar nuclei in the halo. We can 
decompose the observed proton spectrum into its young and old components for a 
given value of K, independent of n^. This has been done by Protheroe (1981) 
for a leakage rate out of the arms which is consistent with the observed 
boron/carbon ratio. 
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e) Other Models 

Secondary positrons have been considered by Stephens (1981) for 
the case of the nested leaky box model of Cowsik and Wilson (1973). In this 
model, cosmic ray sources are surrounded by dense regions of matter in which 
the cosmic rays are partially trapped before leaking out into an outer volume 
where the Sun is located. Escape from the source region is energy dependent, 
resulting in a variation of secondary to primary ratios with energy, while 
escape from the outer region is independent of energy. The effect of the 
matter surrounding the source is to produce a pathlength distribution which is 
deficient in short pathlengths when compared to an exponential (leaky box 
model) distribution. This results in the observed secondary to primary ratios 
(e.g., borori/carbon) being obtained for a lower mean escape length than for, 
e.g., leaky box models. With this lower mean escape length, the predicted 
flux of positrons will be lower than for models with an exponetial pathlength 
distribution, except at the very highest energies. This was indeed the result 
found by Stephens (1981). Other propagation models with a deficiency of short 
pathlengths, e.g., the 'no near sources model' (Lezniak and Webber; 1979) will 
also result in lower positron fluxes than in the leaky box model. 

PROPAGATION OF POSITRONS 


We have made a new calculation of the production spectrum of secondary 
positrons resulting from nuclear interactions of cosmic rays in the inter- 
stellar medium. For the cosmic ray proton spectrum we have taken the range of 
demodulated spectra from the work of Morfill, Volk and Lee (1976). First, we 
calculated the production rate of using fits to the inclusive cross section 
data on !!■*■ production in pp collisions surveyed by Taylor et al. (1976) and 
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supplemented by low energy data of Blobel et al. (1974) and more recent hiqh 
energy data of Guettler et aU (1976) and Johnson et al. (1978). Nuclear 
Interactions involving He, either in the cosmic rays or in interstellar matter 
(assumed to be 10 percent by number), were taken account of as described by 
Giler, Wdowczyk and Wolfendale (1977) using emulsion data of Andersson, 
Otterlund and Stenlund (1979) to scale from pp to pHe interactions. The 
positron production spectrum was then obtained after a full treatment of pion 
and muon decays taking into account the muon decay asymmetry and positrons 
resulting from kaon production (Orth and Buffington, 1976). The production 
spectra obtained for tt+, and e'*’ are given in Figure 1 where the uncertainty 
at low energies due to uncertainties in the demodulation of the proton 
spectrum and at high energies due to uncertainty in the extrapolation with 
energy of the inclusive cross sections are indicated. The result for 
positrons is compared in Figure 2 with those obtained by previous authors and 
found to be in excellent agreement with that of Orth and Buffington (1976). 

Energy losses by synchrotron radiation, inverse Compton interactions, 
bremsstrahlung and ionization are important in determining the shape of the 
positron energy spectrum for a given production spectrum. For synchrotron 
losses, we adopt an r.m.s. iijagnetic field strength of 6 microgauss, the value 
required to give consistency between the observed cosmic ray electron spectrum 
and the radio-synchrotron emission observed from the Galaxy (Rockstroh and 
Webber, 1978; Webber, Simpson and Cane, 1980). This r.m.s. value, is about 
twice as large as that usually adopted for the mean magnetic field strength. 

In addition to the 2. 70 k microwave background, we consider inverse Compton 
scattering off the far infra-red and optical radiation fields. The radiation 
densities we adopt for these fields are 0.47 eV/cm^ and 0.46 eV/cm^ 
corresponding to the local values in the model of Kniffen and Fichtel (1981) 
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which is based on the infra-red survey by Boiss£ et al. (1981) and the stellar 
distributions of Bachall and Soneira (1980). These values lead to ^ (Synch, + 
Compton) " 2.2x10“!^ CeV/s, For Ionization and brerosatrahlung losses we use 
formulae from Ginzburg and Syrovatskii (1964). 

The flux of positrons in the leaky box model (exponential pathlength 
distribution) is given by: 



where PgiE) is the rate of production of positrons (GeV“ls"^q"M, p is the 
density (g cm~3), <t(E)> is the mean cosmic ray age at energy E, and (dE/dt) 
is the rate of energy loss from synchrotron, inverse Compton, bremsstrahlung 
and ionization. 

In the Peters and Westergaard (1977) closed galaxy model, the positron 
flux is made up of two components. The young component is identical to the 
flux calculated for the leaky box model while the old component is obtained 
from equation (3) with <t(E)>-K». The rate of production of positrons in the 
halo depends on the old component of the proton spectrum. We have calculated 
this production rate (shown in Figure 1) for the old component of the proton 
spectrum obtained by Protheroe (1981) for K = 100. Peters and Westergaard 
(1977) found that this value of K was consistent with the observed secondary/ 
primary ratios and it is also consistent with the high energy anti proton data 
(Protheroe, 1981). Since the rate of energy loss depends on density, the old 
component of positron flux will also depend on the density in the halo. 

For the diffusion models, we have used the Monte Carlo technique 
described by Owens and Jokipii (1977b). Analytic treatments are available for 
specific cases of power law injection spectra (eg. Lerche and Schlickeiser, 
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1980). We have used the Monte Carlo technique as the Injection spectrum of 
positrons (Figure 1) Is not a power law, and because this method facilitates 
treatment of a break In the energy dependence of the diffusion coefficient. 


OBSERVED SPECTRA 


In order to reduce systematic differences between the various experi- 
ments, we shall compare our predictions with the observed eV(e‘*‘+e“) ratio 
rather than with the positron spectrum directly. We must therefore consider 
the total Interstellar electron (e’’’ and e") spectrum In some detail. Direct 
measurements have been made up to several hundred GeV; however, below ~ 10 GeV 
the electron spectrum observed rfl»'ect1y differs considerably from the Inter- 
stellar spectrum because of solar modulation. At low energies then the best 
estimates of the Interstellar spectrum may come from radio observations of the 
galactic synchrotron emission. Tan and Ng (1981a) have however recently 
attempted a demodulation of the direct observations and find a local Inter- 
stellar electron density which Is about a factor of 10 lower at 100 Me'/ than 
the spectrum of Webber, Simpson and Cane (1980). This discrepancy will be 
discussed later. We show In Figure 3 a representative sample of the direct 
observations above ~ 5 GeV together with the Interstellar spectrum at low 
energies Inferred from radio data by Webber, Simpson and Cane (1980). The 
Interstellar spectrum we adopt Is shown as the solid line. 

The e+/(e++e") ratio obtained by dividinq the predicted positron flux by 
the observed total electron flux (figure 3) Is plotted In Figure 4(a) for the 
leaky-box, diffusive halo and dynamical halo models and In Figure 4(b) for the 
closed galaxy model (K»100) for various densities In the halo. The observed 
ratios are also given In these figures for comparison. The differences 
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between the predictions shown in Figure 4(a) are small, and we cannot 
distinguish between these models with existing data* From 1 to 10 (5eV, all 
the predictions except for the closed galaxy model with a high density in the 
halo 0.3 cm“^) are consistent with the observations. 

Below 1 GeV none of the predictions fits the observed ratio, hm'ver 
solar modulation must be Considered before drawing conclusions. Tf tne 
modulation is the same for e'*’ and e“, and the modulation can be approximated 
by the force field solution (Gleeson and Axford, 1968) then the observed 
ratios should be shifted to a higher energy corresponding to the observed 
energy plus the mean energy lost in the heliosphere, increasing the 
discrepancy. However, this simple picture of modulation may not he correct 
(Burger and Tanaka, 1970; ilokipii and Kopriva, 1979). In any case, 
demodulation of the data is unlikely to reduce the discrepancy unless 
positrons are modulated differently from electrons. If the cosmic ray 
electron density varies over distances of ~ 100 pc, then the local 
interstellar electron spectrum may be lower than that obtained from radio 
data. This has been suggested by Strong and Wolfendale (1978) and Tan and Mg 
(1981b) and may account for the discrepancy. Alternatively, the mean escape 
length of electrons or positrons may differ from that of nuclei (Giler, 
Wdowczyk and Wolfendale, 1977). 

The models discussed earlier showing a deficiency of short pathlengths in 
the pathlength distribution, e.g., the nested leaky box model, give a lower 
positron flux in this energy range and hence give a worse fit to the data. 
Motivated by the high energy antiproton data (Golden et al., 1979), Cowsik and 
Gatsser (1981) have, however, suggested a modification to the nested leaky box 
model which can give enhanced antiproton production in the galaxy without 
affecting the secondary/primary ratios. This modification, which could be 


12 


applied to most of the other propagation models as well, Involves the addition 
of a set of cosmic ray sources shrouded with ~ so g/cm^ of matter. Cosmic ray 
nuclei would interact on traversing the matter, producing pions, antiprotons, 
etc., and spall, eliminating the complex nuclei (i.e., those heavier than 
protons). The neutral pions would decay into gamma-rays j these additional 
sources are thus to be identified with the discrete galactic gamma-ray 
sources (Swanenburg et al., 1980). Positrons would result from the positive 
pions which are produced and may or may not contribute importantly to the 
cosmic ray positron flux depending on the strength of the magnetic fields 
associated with these sources. In any case, the energy spectrum of these 
additional positrons as seen at Earth would be steeper than for those produced 
in the interstellar medium because of energy losses both in the source region 
and on traversing the finite distances from the sources to the Earth. The 
addition of such a component may possibly improve the agreement between the 
predicted and observed fluxes at low energies. 

Above 10 GeV the observed ratio lies above the predictions except for the 
closed galaxy models having a high matter density in the halo. The statis- 
tical errors for these data are however large, and new measurements are 
required before conclusions can be drawn about possible primary origin. 


CONCLUSIONS 


The production spectrum of secondary positrons has been calculated over 
the energy range 100 MeV to I TeV. Observations of the positron spectrum in 
this energy range should provide information about the propagation of cosmic 
rays in the galaxy and solar modulation. In particular, they may enable us to 
distinguish between the various propagation models that have been proposed. 
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With f;he present measurements of the eV(e‘*' + e") ratio, we are unable to 
distinguish between leaky box, diffusive halo and dynamical halo models. For 
progress at energies below a few fieV, a greater understanding of the solar 
modulation of electrons and positrons and the relationship between the l ocal 
interstellar spectrum and the observed radio data is required. In addition, 
new experiments with higher exposure factors will be required as well as 
improved measurements of the boron/carbon ratio and l^Re abundance to 
constrain the propagation models. 

The data rule out a large primary positron component at high energies 
distributed uniformly throughout the Galaxy. A component as large as ~ 2 
percent of the observed electron spectrum is, however, allowed within the 
present uncertainties. The observation of a gamma-ray line at 0,511 MeV * 
(Leventhal, MacCullum and Stang, 1970) has been interpreted to indicate that 
low energy positrons are copiously produced in the galaxy (Ramaty and 
Lingenfelter, 1979). A primary positron component as large as a few percent 
could arise if only a small fraction of these were accelerated to high 
energies (Lingenfelter and Ramaty, 1979). Definitive statements about primary 
positrons must however await new measurements. 
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Figure Captions 


Figure 1; Production rates of and e'*' per interstellar nucleon in the 

disk of the gala)^/? Uncertainties associated with the demodulation of the 
cosmic ray proton spectrum and with extrapolation of cross sections to high 
energies are indicated by hatching. In addition, the production spectra are 
uncertain by at most a further 15 percent due to uncertainties in the 
tranverse momentum distribution of pions produced in pp interactions. Also 
shown, e'*' (old component), is the production rate of positrons in the halo of 
the closed galaxy model for K=100. 

Figure 2; Production rate of e+ per gram of interstellar matter from the 
present work is compared with previous results. For other references to 
earlier work, see Orth and Buffington (1976). 

Figure 3; A representative sample of electron spectrum measurements. The 
cosmic ray electron spectrum used in the present work is indicated (solid 
line). 

Figure 4; Comparison of observed e+/(e^+e") ratio with those obtained by 
dividing predicted e"** flux by observed (e^ + e“) flux for: (a) leaky box, 
diffusive halo and dynamical halo models; (b) closed galaxy model (K=inn) for 
various densities of neutral matter in the halo (dashed line for ionized 
matter). Data are from: Buffington, Orth and Smoot (1975) (■); Daugherty, 
Hartman and Schmidt (1975) (•); Fanselow, et al. (1969) (D); Hartman and 
Pellerin (1976) (O). The error bar attached to the prediction for the 
diffusive halo model indicates the precision of our Monte Carlo calculations. 
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